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Abstract 

We show how nonzero Majorana neutrino masses can be radiatively generated 
by extending the MSSM with leptoquark chiral multiplets without violating 
R-parity. It is found that, with these particles, the R-parity conservation 
does not imply lepton number conservation. Neutrino masses generated at a 
one-loop level are closely related to the down quark mass matrix. The ratio 
of neutrino mass-squared splittings Am^^/Am^ obtained is naturally 
close to Am^_ d / Aml_ s ~ 10~ 3 which is in the right region required to explain 
both the atmospheric neutrino data and the MSW solutions for the solar 
neutrino data. 
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I. INTRODUCTION 



Many recent experimental observations have hinted at non-vanishing neutrino masses 
[IJ. The strongest hint is from the Super-Kamiokande underground experiment 0. They 
observed an up-down asymmetry of atmospheric neutrino flux, suggesting an oscillation 
of into v T (or u s ). The mass squared splitting is found in the range of 1.5 x 10~ 3 < 
Am^.^ < 6 x 10 _3 eV 2 (with 95% C.L.), and the mixing angle is close to the maximal 
mixing H. 

In addition, the observed solar neutrino fluxes [Q are well below the Standard Solar Model 
predictions. The details may be explained by invoking neutrino mixing between v e and v^. 
There are two types of solutions. One is the solution due to the MSW mechanism || with 
either large or small mixing angles ||. For the small mixing angle (SMA) solution, one needs 
Aml 2 _ Ul = (0.4-1) x 10~ 5 eV 2 , and sin 2 29 = (0.12-1.2) x 10~ 2 . For the large mixing angle 
(LMA) solution, the mass splitting and mixing angle are Am 2 2 _^ = (0.6 — 20) x 10~ 5 eV 2 , 
and sin 2 29 ~ 0.76. The other one is the so-called just-so solution, which corresponds to 
a neutrino oscillation in vacuum [[?[]. For this case, one needs the mass-squared splitting 
Am 2 = 8.0 x 10~ 11 eV 2 and the mixing angle sin 2 29 = 0.75 for the best fit solution ||. 

If neutrinos have masses, they are expected to be small from laboratory bounds and 
cosmological consideration |8"|-|TT]|] . In the Standard Model (SM) neutrinos are massless. 
One needs to extend the SM to accommodate massive neutrinos. A neutrino mass can be 
generated by a Higgs Yukawa interaction like any other fermions if right-handed neutrinos 
are introduced. However, this corresponds to an unnaturally small Higgs Yukawa coupling, 
A ~ m u /v < eV/246GeV ~ 10~ n , where v is the vacuum expectation value of the Higgs 
field. There have been extensive studies of neutrino masses in models beyond the SM 



n|]12|| . To explain the smallness of the neutrino mass, a popular approach is to use the see- 
saw mechanism to trade a light neutrino with a super heavy neutrino |13| or to radiatively 
generate neutrino masses such that the smallness of neutrino masses are due to the loop 



suppression 
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In studying neutrino masses, it is interesting to note that the ratio of the mass- squared 
splittings Aml 2 _ Ui /Aml 3 _ U2 for the MSW solutions to the solar neutrino problem can be 
close to Aml_ d / Aml_ s ~ 10~ 3 . This may be regarded as a hint that the neutrino mass 
matrix is proportional to the down quark mass matrix. In this case the heaviest neutrino 
has a mass less than ~ 0.08 eV. It is interesting to see if such a mass matrix can be obtained 
in a natural way. In this paper we study a model in which the Majorana neutrino mass is 
radiatively generated by supersymmetric leptoquarks and the mass matrix is closely related 
to the down quark mass matrix. 

The supersymmetry (SUSY) is one of the leading candidates for physics beyond the SM 
15fl . In order to generate neutrino masses one needs to violate R-parity and/or introduce 



new physics [fL6|,|17 |. We will show that it is possible to generate Majorana neutrino masses 
without the violation of R-parity if leptoquarks are introduced in the model. In the presence 
of these leptoquarks, the lepton number (L) can be broken explicitly even if the R-parity is 
conserved. The lepton number violating interaction can generate Majorana neutrino masses 
at a one-loop level and naturally lead to a mass matrix proportional to the down quark mass 
matrix. This is due to a special feature of the scalar leptoquark interaction. It can be shown 
that a lepton number violating vector leptoquark model will relate a neutrino mass matrix 
to the up quark mass matrix, which would not naturally give neutrino mass hierarchies 
mentioned above. In the SUSY extension with the SM gauge group, vector leptoquarks are 
not allowed. This provides additional motivation to study SUSY models with leptoquark 
multiplets. 



II. THE MODEL 

The matter and Higgs superfields in the MSSM are L, E, Q, U, D, Hi and H2. Their 
quantum numbers under SU(3) c xSU(2) L xU(l)y are (1, 2, -1/2), (1, 1, +1), (3, 2, +1/6), 
(3*, 1, -2/3), (3*, 1, +1/3), (1, 2, -1/2) and (1, 2, +1/2), respectively. To avoid lepton 
number and baryon number violations, R-parity is introduced. The ordinary fields, the 
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quarks, leptons and Higgses, are assigned with positive R-parity and their super-partners 
are assigned with negative R-parity. There are in general seven types of leptoquarks which 
can couple directly to the quark and lepton field products, U x E, Q x L, D x E, U x L, 
Q x E, and D x L, at the tree level with their gauge quantum numbers given, respectively, 
by 

S : (3, 1, -1/3), S' : (3% 1, +1/3), f : (3*, 3, +1/3), S" : (3, 1, -4/3), 

R : (3, 2, +7/6), R> : (3*, 2, -7/6), I : (3, 2, +1/6). (1) 

The R-parity assignments of these leptoquark supermultiplets are determined by their in- 
teraction with lepton and quark multiplets to be the same as that of Higgs multiplets. 

For the purpose of generating Majorana neutrino masses only two types of leptoquarks 
are needed which are S' and /. These leptoquarks provide two units of the lepton number 
violation needed for Majorana neutrino masses through the R-parity conserving interaction 
in the superpotential by a term, eijH\VS' . In order to cancel the gauge anomalies due to 
the fermion partners of leptoquarks for consistency more leptoquark multiplets are required. 
We note that the quantum number of S is conjugate to that of S' and accordingly it can be 
used to cancel the anomalies due to S'. We need to introduce a new leptoquark multiplet 
I' with quantum number conjugates to J to cancel the anomalies. We also note that the 
R-parity does not forbid a baryon number violating term UDS'. To avoid rapid proton 
decay, we must either assume the coupling is very small, or introduce further symmetry. To 
make use of the symmetry, we may assign a negative Z2 parity to U, D and Q, while all 
other matter fields in the MSSM have a positive Z2 parity. The Z2 parity of leptoquarks 
is negative as determined from leptoquark interactions. Thus, the Z 2 symmetry will forbid 
the baryon number violating terms. In short, we are considering the leptoquark interaction 
instead of the diquark interaction. 
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Superfield 


SU(3)xSU(2)xU(l)xL 


Boson Fields 


Fermionic Partners 


S 


(3 1 -1/3- -1) 


St 


Si 


S' 


(3*, 1, +1/3; +1) 




Si 


i 


(3, 2, +1/6; +1) 


(Ol,Pl) 


(Ol,Pl) 


P 


(3*, 2, -1/6; -1) 


{P*R, -Or) 


{pi -oi) 



Table 1. The particle content of leptoquark chiral supermultiplets related to the neutrino 
mass generation. The L denotes the lepton number. 

The model we now have is an extended MSSM with the standard particle content plus 
leptoquarks and leptoquarkinos in Table 1. With these multiplets we can construct the 
leptoquark R-parity conserving superpotential Wlq with 

W LQ = -(h*) ab U a E b S - {h') ab e l3 Q\L{S' - (h") ab e t ,D a LlP 

+g e t] H\PS' + M s S d S' d + Mje tJ P l P, (2) 

where a and b are generation indices of quarks and leptons in the weak basis. Note that Sr 
and Pl can be mixed by the g term when the Hi takes its vacuum expectation value (VEV) 
V\\ otherwise, the masses are Ms and Mj respectively. 

With this superpotential and the Higgs VEVs, we get leptoquark mixings of Sl with P L , 
Sr with Pr, and Sr with Pl. The F-term causes mixings, and the D-term contributes to 
diagonal parts of the mass matrix only. The mixing terms are 

V = gv x (M s S{P L + Mj S R P R -fitm/3 S R P L ) + h.c. (3) 

There are more terms due to the SUSY soft breaking, 




+M 2 Sr _ Sl S r S l + Ml r e l3 n>J + Mmse^HlPSi + h. c. (4) 

Note that the (standard) SUSY soft breaking terms do not include H\IS L or H\lS R and 
so on [|T5| . The mass parameters Mj in the Eq. (^) are naturally in the SUSY breaking 
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scale. The MHis^ijH\P S R term violates lepton number by two units. After spontaneous 
symmetry breaking this term produces the mixing of Sr and which plays a crucial 
role in generating Majorana masses at a one-loop level. This is the dominant source for the 
lepton number violation related to the Majorana neutrino masses although other leptoquarks 
contribute through mixings among Sl, Sr, Pl and Pr. We note that the mass parameters, 
Ms and Mj, from the superpotential are free parameters which may in general be smaller 
or larger than the SUSY breaking scale Msusy, while those mass parameters contained in 
soft breaking terms are naturally in the scale of Msusy- 



III. RADIATIVE MAJORANA NEUTRINO MASSES 

The Majorana neutrino mass term 1/2 {v Ra m a b v^ b + h.c.) can be obtained from one-loop 
diagrams with all possible leptoquark mass insertions in the internal leptoquark line and the 
internal quark line. The leading contribution comes from a mass insertion each in the quark 
line and the other insertion in the leptoquark line as shown in Figure 1. The mass matrix 
m a b is symmetric and the leading term is given by, 

m ab = ^^^M D h K + h" T M T D h'U (5) 

where N c = 3 is the number of color, is the down type quark mass matrix, Mlq is 
the average leptoquark mass and AJU 2 = -Mp s = — g^i/itan/? + MrisVi ~ Msusy v 
denotes the Pl — Sr mixing element as shown in Eq. (3) and Eq. (4). 

The quark mass matrix can be diagonalized in the usual way, VdiM D Vj R = M^ iag , where 
the diagonal mass matrix is given by [[HJ . 



^D iag (/U = diag(m d , m s , m b ) 

= m b (ii) diag((1.36 ± 0.42) x 1(T 3 , (2.66 ± 0.87) x 10~ 2 , 1 ± 0.05) . (6) 



From Eq. (|^), it is clear that the neutrino mass matrix obtained is closely related to the down 
quark mass matrix. We will work in a quark field basis where the down quark mass matrix 
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is hermitian, which can always be done |21| . If small CP violation entries are neglected the 
down quark mass matrix is symmetric, V^l becomes orthogonal and equals to V^r- This 
choice of basis is convenient for our purpose since the Majorana neutrino mass matrix is 
symmetric. 

To finally obtain numerical numbers for neutrino masses, we need to know parameters 
related to leptoquarks. The couplings h', h" and the leptoquark mass Mlq relevant to 
neutrino masses in our model are unknown which must be subjected to experimental con- 
straints. The typical experimental lower bound on the leptoquark mass is about 200 GeV 



22fl . There are many studies on leptoquark couplings and masses from indirect searches 
23|-p5|. We now consider the constraints on h' ab , h!' ah and M L q. We follow ref. ( p3| ) to 



evaluate te constraints. We note that h', h" are denoted as \ls d i 2 m re f- [0 



For h' ab the strongest constraint comes from the lepton flavor- changing rate for a muon 
to turn into an electron when scattered off a nucleus, such as /xTi — > eTi. This gives the 



constraint on the coupling of Sr as |23|| 



^1^2/^LQ^O.OlTeV- 1 . (7) 



Similar consideration gives a constraint on the coupling of P with y h'dh'^/MLQ < 
0.01 TeV -1 . The leptoquark exchanges also contributes to the kaon semileptonic decay, 
K + — > 7i + uu. By requiring that the prediction is smaller then the observed bound of 



Br(K+ -> 7r + z/^) = 0.42±§;§£ x 10" 9 gives constraints ||||2§, 



h' xl h' 2V /M LQ , ^KKvIMlq < 0.03 TeV- 1 . (8) 

More stringent constraints on the P couplings come from lepton family number violating 
leptonic kaon decay, K° L — > /i ± e =F . Using the constraint in ref. p3| and the new experimental 
result 0, Br{K° L -> ^e T ) < 3.3 x 10 -11 , we have 



K^/Mlq, ^h>{ 2 h^/M LQ < 6 x 10^ 3 TeV" 1 . (9) 

One can find constraints on other elements of h' ab , h!' ah in ref. |]23|| , they are in general weaker 
than the ones discussed here. 



As we turn on the mixing of Sr with Pl, these chiral leptoquarks become non-chiral 
leptoquarks, there are further constraints for them. The constraints from the chirally sup- 
pressed pion decay, 7r + — > e + ui give |8|| 



h'nK,^^, < 4 x 10-TeV-i. (10) 

LQ LQ 

The similar kaon leptonic decay gives a weaker constraint due to m s > m u , m&. The mixing 
of Sr with Pl contributes to the lepton number violating (AL = 2) K + — > -n^vivv decay. 
This gives 



K,h'i,^^, < 2 x 10-TeV-. (11) 

LQ LQ 

The factor y/ M S usyv / M L q in Eqs. ([U]) and ([11]) reminds us that these non-chiral leptoquark 
behavior is due to the mixing and the constraints are weaker than the constraints obtained 



in Eq. <H) if V Msusyv / Mlq << 1. Without any mixing, there is no non-chiral leptoquark 
in the SUSY model. 

To estimate what value of m a b may be, we may work with a simplified situation for 
illustration, with h'^ ~ h"j ~ h5ij. We obtain a very simple and very interesting mass 
matrix with the mass eigenvalues of the neutrinos given by 

N c AM 2 , 2/ v , , 

[m vi , m V2 , m U3 ) = - h [m d , m s , m b ), (12) 

and the mixing matrix is given by VdL- We note that the small mixing behavior of Cabbibo- 
Kobayashi-Maskawa matrix (Vckm — V u lV^ l ) does not necessarily imply any small mixing 
pattern in the down quark mass matrix or up quark mass matrix separately. There are 
parameter spaces in which the required mixing angles for the atmospherical and the solar 
neutrino data can be produced. 

To obtain AmJ rJ/2 to be below 6 x 10~ 3 eV 2 as required by the atmospheric neutrino 
data, H/Mlq is restricted to be less than 0.55 x 10 _4 TeV _1 which is safely below the above 
constraints for Mlq > Msusy- We note that for the particular case where hi and h" are 
universal and flavor blind, constraints in Eqs. (H), @, fllCf ) and (|iT|) all apply. Using 

8 



h/MiQ = 0.53 x 10 4 TeV 1 with m& at the top quark mass scale and Eq. fl5|), we obtain 
Majorana neutrino masses, 

m Vl = (9.92 ±3.08) x 10~ 5 eV, 
m U2 = (1.94 ±0.64) x 10~ 3 eV, 

m U3 = (7.30 ± 0.36) x 10~ 2 eV. (13) 
These values correspond to mass-squared splittings of 

Aml 3 „ U2 = (5.32 ± 0.53) x 10" 3 eV 2 , 

Aml 2 _ Ui = (3.75 ± 2.48) x 10~ 6 eV 2 , (14) 

which are in the right regions with the atmospheric neutrino solution and the MSW solar 
neutrino small mixing angle and large mixing angle solutions. These Majorana neutrino 
masses are well below the bound on the effective Majorana mass (m.„ e „ e < 0.2 eV) from the 
negative result of the neutrinoless double beta decay |J. 

If we use vector leptoquarks instead of scalar leptoquarks, we would have obtained neu- 
trino masses proportional to up type quark masses. There is an important difference when 
compared with the scalar leptoquark case discussed above. One would naturally obtain 
) to be proportional to (m u , m c , m t ). The ratio Am^^/Am^.^ is then 
given by Am 2 _ c /Am 2 _ i ~2x 10~ 5 which disagrees with the data. 

IV. CONCLUSION 

In this paper, we have shown that a simple extension of the MSSM with leptoquark 
multiplets contains a lepton number violating interaction which can generate Majorana 
neutrino masses at a one-loop level without the breaking of R-parity. One may contrast this 
scenario to the R-parity violating extension of MSSM. In both cases, in order to obtain a 
light neutrino mass, one needs the ratio h/M (or X^/Mmssm in the extended MSSM) to be 
of the order of 10~ 4 TeV™ 1 . In the particle content of the MSSM all particle masses (M MSS m) 
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should be below Msusy, and this would correspond to suppressed R-parity violating Yukawa 
couplings (Ag). In our case, in the presence of the leptoquark multiplets, the Yukawa 
couplings can in general be not artificially suppressed. The AL = 2 interaction in this 
model is favored than the AL = 1 interaction in the extended MSSM in light of leptogenesis 
FT] ]. Furthermore, the radiatively generated neutrino mass matrix in this model is closely 



related to the down quark mass matrix as leptoquarks couple to leptons and down quarks. In 
the simplest scenario with the coupling matrices for leptoquarks, quarks and leptons being 
proportional to unit matrix, that is, h'^ ~ h"j ~ hSij, it predicts a hierarchy structure in 
neutrino masses Am^.^/Am^^ ~ Am 2 d _ s / Am 2 s _ b which provides mass-squared splittings 
required to explain both the atmospheric neutrino and the MSW solar neutrino data. 

More complicated structures for h' and h" can lead to different mass hierarchy patterns. 
The simple neutrino mass matrix that is proportional to the down quark mass matrix, 
depends on that the matrices h'(h") being proportional to the unit matrix. It remains to be 
seen that if this can be obtained from some symmetry principle. 

This work was supported in part by a grant from National Science Council of Republic 
of China grants NSC88-2112-M002-001Y and NSC88-2112-M002-041. 
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FIGURES 




FIG. 1. The leading one-loop diagram which gives rise to nonzero neutrino masses, the X 



represent a mass insertion. 
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